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Abstract 

The evolution equations of the vorticities of the electrons, ions and photons in a pre- 
decoupling plasma are derived, in a fully inhomogeneous geometry, by combining the general 
relativistic gradient expansion and the drift approximation within the Adler-Misner-Deser 
decomposition. The vorticity transfer between the different species is discussed in this novel 
framework and a set of general conservation laws, connecting the vorticities of the three- 
component plasma with the magnetic field intensity, is derived. After demonstrating that 
a source of large-scale vorticity resides in the spatial gradients of the geometry and of the 
electromagnetic sources, the total vorticity is estimated to lowest order in the spatial gra- 
dients and by enforcing the validity of the momentum constraint. By acknowledging the 
current bounds on the tensor to scalar ratio in the (minimal) tensor extension of the ACDM 
paradigm the maximal comoving magnetic field induced by the total vorticity turns out 
to be, at most, of the order of 10 -37 G over the typical comoving scales ranging between 
1 and 10 Mpc. While the obtained results seem to be irrelevant for seeding a reasonable 
galactic dynamo action, they demonstrate how the proposed fully inhomogeneous treatment 
can be used for the systematic scrutiny of pre-decoupling plasmas beyond the conventional 
perturbative expansions. 
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1 Primordial vorticities? 



The observed Universe might originate from a strongly coupled electromagnetic plasma exist- 
ing prior to photon decoupling where the angular momentum transfer between ions, electrons 
and photons in an expanding space-time geometry leads to the formation of large-scale vor- 
tices as speculated by various authors including, with slightly different perspectives, Hoyle 
PQ , Harrison [2j [3j 0] , Mishustin and Ruzmaikin [5] , Ozernoy and Chernin [61 [7J |8] and oth- 
ers (see also |9j). The primordial vorticity, if present in the pre- decoupling plasma, might 
lead eventually to the formation of large-scale magnetic fields possibly relevant for galactic 
magnetogenesis. 

The physical description of the angular momentum exchange between ions, electrons and 
photons can be realized by appropriately translating the evolution equations describing ion- 
ized gases [TQl E] to an expanding geometry supplemented by its own relativistic fluctuations 
[T2"|IT5]. In the context of the ACDM paradigrro, it is both reasonable and justified to assume 
that the background geometry is conformally flat and that its inhomogeneities stem from the 
relativistic fluctuations of the spatial curvature described either in gauge-invariant terms or 
in an appropriate gauge. The latter assumption rests exactly on the absence of large-scale 
vorticity which is assumed to be vanishing at least within the current observational precision. 

A gross argument could suggest that the vorticity must be negligible for ACDM initial 
conditions, since it is the curl of a velocity. Thanks to the momentum constraint (connecting 
the first derivatives of the linearized fluctuations of the geometry to the peculiar velocities), 
the total velocity field is subleading when compared with the density contrasts or with the 
curvature perturbation for typical scales larger than the Hubble radius and in the case of the 
conventional adiabatic initial conditions postulated in the vanilla ACDM scenario. The latter 
argument suggests that the treatment of large-scale vorticity assumes, more or less tacitly, 
a correct treatment of the spatial gradients. To transform this incomplete observation in a 
more rigorous approach it is necessary to introduce a description of the vorticity which does 
not rely on the purported smallness or largeness of the gravitational fluctuations. It is rather 
desirable to describe the angular momentum exchange between ions, electrons and photons in 
a gravitating plasma which is also fully inhomogeneous. By fully inhomogeneous plasma we 
mean the situation where not only the concentrations of charged and neutral species depend, 
in an arbitrary manner, upon the spatial coordinates but where the geometry as well as the 
electromagnetic fields are not homogeneous. It has been recently argued [H] that such a 
description can be rather effective for the analysis of a wide range of phenomena including 
the physics of pre-decoupling plasmas. In the present paper the results of Ref. [H] shall 
be first extended and then applied to a concrete situation with the purpose of obtaining an 
explicit set of equations describing the evolution of the vorticities of the various species of the 
plasma. The proposal of [2] is built on the fully inhomogeneous description of the geometry 

3 The acronym ACDM (where A denotes the dark energy component and CDM stands for cold dark 
matter) and the terminology concordance paradigm will be used interchangeably. 



2 



in terms of the Adler-Misner-Deser (ADM) variables [T5j ITB] which are customarily exploited 
for the implementation of the general relativistic gradient expansion [T7J [HJ [191 1201 IZD 122] . 
The second key ingredient of Ref. [H] is the fully inhomogeneous description of cold plasmas 
in flat space which is the starting points of the analysis of nonlinear effects in kinetic theory 
and in magnetohydrodynamics (see, e.g. [231 [2U 125] ) • Consequently, the vorticity exchange 
between ions, electrons and photons can be analyzed in gravitating plasmas with the help 
of an expansion scheme which involves not only the gradients of the geometry, but also the 
gradients of the electromagnetic sources, by so combining the general relativistic gradient 
expansion and the drift approximation (sometimes dubbed giuding center approximation) 
typical of cold plasmas. 

The approach pursued in this paper reproduces, in the conformally flat limit, the conven- 
tional treatment which will be made more precise in section [2j The evolution of the vorticity 
in gravitating plasmas which are also fully inhomogeneous will be discussed in section [3j In 
sections H] and [5] the total vorticity of the geometry will be computed within the gradient 
expansion and estimated in the framework of the ACDM paradigm. The maximal magnetic 
field induced by the total vorticity will be computed in section [61 Section [7] contains our 
concluding remarks. In appendix [A] some useful complements have been included to make 
the paper self-contained while in appendix IB1 useful details on the calculations of correlation 
functions of multiple fields in real space have been included for the technical benefit of the 
interested readers. 

2 Vorticities in conventional perturbative expansions 

The treatment proposed here differs slightly from the one of Refs. [21 El SI E] for three 
reasons: (i) the conformal time coordinate is preferred to the cosmic time; (ii) the relativistic 
fluctuations of the geometry are included in the longitudinal gauge; (iii) the three-fluid, 
two-fluids and one fluid descriptions are discussed more explicitly within the appropriate 
temperature ranges where they are applicable. 

The conformal flatness of the geometry does not imply the invariance of the system under 
the Weyl rescaling of the metric. Such a potential symmetry is broken by the masses of the 
electrons and ions which are crucial in the large-scale evolution of the vorticity. The con- 
siderations of the present section can also be formulated in the case of a geometry which is 
not spatially flat; this is not essential since the subsequent generalizations will automatically 
include also geometries which are not necessarily spatially flat. Consider first the case of a 
conformally flat background geometry characterized by a metric tensor = a 2 (r)?7 At j y and 
supplemented by the corresponding relativistic fluctuations which we write in the longitudi- 
nal gauge 

5 s 9oo{x, t) = 2 a 2 (r)(j)(x, r), 6 B gij(x, r) = 2 a 2 (r)^(x, r)<% (2.1) 
note that 5 S describes a metric perturbation which preserves the scalar nature of the fluctu- 
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ation since, in the ACDM paradigm, the dominant source of inhomogeneity comes from the 
scalar modes of the geometry. By defining the comoving electromagnetic fields E and B as 
well as the comoving concentrations of electrons and ions (i.e. n e and n{) 

E(x, t) = a 2 {r)S{x, r), B(x, r) = a 2 (r)B{x , r), 

7ii(x, t) = a 3 (r)hi(x, r), n e (x,r) = a 3 (T)h e (x,r), (2.2) 

Maxwell's equations read 

V- J E = 47re(n i -n e ), V • B = 0, (2.3) 
Vx£ = -<9 r £?, Vx5 = 4vre(ni - n e tf e ) + d T E. (2.4) 

In Eqs. (12. 2p . (12.31) and (12.41) all the fields are appropriately rescaled so that the resulting 
equations are formally equivalent to the ones of the flat space-time. The peculiar velocities 
of the ions, electrons and photons obey the following set of equations^ 



CTl -* -» -* 4 

d T v e + Hv e = e -[E + v e x B] - V0 + -^aT JC (v y - v c ) + aT^ - v c ), (2.5) 

Pe A 4 3 p e 

T Vi + Hv- l = ^j[E + ^ x B] - V0 + ^ar 7i (t7 7 - v,) + aT ci ^(v e - v{), (2.6) 
Pi a 4 3 pi pi 

d T Vj = — -V<5 7 — V0 + ar 7 i(-Ui — Uy) + aT je (v e — (2.7) 
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In Eqs. (I2.5p -( 12T7|) the relativistic fluctuations of the geometry are included from the very 
beginning in terms of the longitudinal gauge variables of Eq. (12. ip ; the electron-photon, 
electron-ion and ion-photon rates of momentum exchange appearing in Eqs. (I2.5p - (l2.7p are 
given bj@: 



e 2 \ 2 8 fe 2 ^ 2 



r 7e = n e a ei , r 7i = nicr i7 , cx e7 = -7r — , cr i7 = -n[ — , (2.8) 

3 \m e J 3 \mj 

fT~ e 4 3 / T 3 

T ei = n c J — a ci = r ie , o- ci = — lnA c , A c = — J- — . (2.9) 
V m e 1 z 2e d V ^-e 71 " 

Note that, in Eq. (12. 8p and (12. 9p . T and n are, respectively, physical temperatures and 
physical concentrations. If the rates and the cross sections would be consistently expressed 
in terms of comoving temperatures T = aT and comoving concentrations n = a 3 n the 
corr esponding rates will inherit a scale factor for each mass. For instance aT e i becomes 
n c Jt /(m c a) (e 4 /T 2 ) lnA c , if comoving temperature and concentrations are used. 

Let us then define the vorticities associated with the peculiar velocities of the various 
species 

u e (x, t) = V x v e , uji(x,r) —Vxvi, uj 7 (x,t) = V x t7 7 , (2-10) 



4 As usual H = d T In a and its relation with the Hubble rate is simply H = aH. 
5 Note that T denotes the temperature and Ac is the Coulomb logarithm [Ttil 111] , 
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and their corresponding three- divergences: 

9 c (x, r) = V • v e , t) = V • Vi, 9j(x, r) = V • u 7 . 



(2.11) 



The evolution equations of the vorticities and of the divergences can be obtained by taking, 
respectively, the curl and the divergence of Eqs. (I2.5p - (j2.7p and by using Eqs. (12. 3p and 
(12. 4p . To simplify the obtained expressions it is useful to introduce the total comoving charge 
density and the comoving current density. 

p q = e(n i -n e ), J — efav-y — n e v e ). (2.12) 

Thus, the evolution of the vorticities and of the divergences of the electrons are, respectively, 



d T co c + T-Luj c 



d T 6 e + H6 C 



p e a q 
4p 7 



d T B + {v c ■ V)B + 9 C B - (B ■ V)v c 



+ -—aT ie (uj - u e ) + aT ei (ui- u e 

3 Pe 



+ 



p c a 
4^ 

3 Pe 



1 

aT 



47rp q + u e ■ B — A-Kv c ■ J — v c ■ d T E 

7 e(^7 — #e) + oTei(0i ~ #e), 



- V" 



(2.13) 



(2.14) 



Conversely the vorticity and the three-divergence of the ions evolve as: 



d T Ui + H tOi 



d T 6i + H8i 



em 



Pi or 



d T B + (v { ■ V)B + 9iB - {B ■ V)^ 



+ -— ar 7i (a; 7 - Uj) + aT ci — (u c - £{) 
3 Pi Pi 



eni 
Pi a 4 

3 Pi 



47rp q + uj\ ■ B — 4iiVi ■ J — v- x ■ d T E 
P 



-V 2 



+ -— ar 7 ;(6' 7 — 9\) + aT ei 



Pi 



-ex 



(2.15) 



(2.16) 



Finally, the evolution equations for the photons are given by: 



d T Cu 1 
d T 6^ 



-^v 2 5 7 - v 2 ^ + ar 7i (#i - e 7 ) + ar 7e (# e 



(2.17) 
(2.18) 



ar 7 ;(a;i — w 7 ) + aT ie (u e — w 7 ), 
1 
4 

The system described by the set of equations deduced so far will be considered as glob 
ally neutral. In particular, prior to photon decoupling, the electron and ion (comoving 
concentrations have a common value no, i.e. n\ = n e = no wher^ 



no = Vbon^, 



r/ b o = 6.177 x 10 



.10/ \ /2.725K\ 3 



0.02258 



T. 



7 o 



(2.19) 



6 If not otherwise stated the pivotal values of the cosmological parameters will be the ones determined 
from the WMAP 7yr data alone in the light of the ACDM paradigm. 
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and T 7 q is the present value of the CMB temperature determining the concentration of the 
photons; fibo is the present value of the critical fraction of baryons, while ho is the Hubble 
constant in units of 100Km/(Mpc x sec). The system of Eqs. (12.13]) — (j2.18p is coupled with 
the evolution of the density contrasts of the electrons, ions and photons (i.e. 8 e , 5i and <5 7 ) 

S' e = -9 C + 3^' \E ■ v e , 5{ = -9, + 3^' + -\E ■ v u (2.20) 

5; = Aif,' - U T (2.21) 

Finally the metric fluctuations, the density contrasts and the divergences of the peculiar 
velocities are both determined and constrained by the perturbed Einstein equations (see, 
e.g. Eqs. (2.43)-(2.46) in the first article of Ref. [12J). Concerning the system of Eqs. 
(I2.13p -f l2.18p two comments are in order: 

• Eqs. (12"33]) - (I2TT4]) (as well as Eqs. (l2TT5]) - (l2Tll)]) ) couple together the evolution of the 
vorticities, the evolution of the divergences and the gradients of the magnetic field; 
while in the linearized approximation the spatial gradients are simply neglected, in the 
forthcoming sections the evolution of the vorticity will be studied to a given order in 
the spatial gradients; 

• the electron and ion masses break the Weyl rescaling of the whole system of equations; 
this aspect can be appreciated by noticing that the prefactor appearing in front of the 
square brackets at the right hand side of Eqs. (]2TT3]) - (f2TTI]) and Eqs. (12TT5) -( 12TT5]) is, 
respectively, e/(m e a) and e/{mia). 

Equations fl2.13p - fl2.18l) have three different scales of vorticity exchange: the photon-ion, 
the photon-electron and the electron ion rates whose respective magnitude determines the 
subleading terms and the different dynamical regimes. By taking the ratios of the two rates 
appearing at the right hand side of Eqs. (12 . 13[) and (12.151) the following two dimensionless 
ratios can be constructed^: 

3p e r ei 135C(3) (Ty^ /T\~ 5 / 2 

~a — f = ifi 5 — r/ bo lnA c = — , (2.22) 

4p 7 l 7e l0 7T° \m e J V-ter/ 

3p c r ei (m p \ 2 ( T \-V 2 _ / T \- 5 / 2 



4 p 7 T 7 i \TTl c J \X^ 7 / V^~i 



(2.23) 



where ((3) = 1.202... and the ion mass has been estimated through the proton mass; the 
effective temperatures T e7 and T i7 introduced in the second equality of Eqs. ( 12.221) and ( I2.23j) 
are defined as: 

T e7 = meA/" 2 ^ 5 , T i7 = m-^V^ N = ^^lnA c . (2.24) 



r Note that p x must simplify when taking the ratio of the two rates in Eq. (|2.15[) . 
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In explicit terms and for the fiducial set of cosmological parameters determined on the basis 
of the WMAP 7yr data alone in the light of the ACDM scenario 



T e7 = 88.6 f -^f) eV, T h = 36.08 ( _M^MlV /5 keV . (2.25) 
7 V0.02258/ ' 7 V0.02258/ v ; 

On the basis of Eq. (I2.25j) there are three different dynamical regimes. When T > T i7 the 
ion-photon and the electron-photon rates dominate against the Coulomb rate: in this regime 
the photons, electrons and ions are all coupled together and form a unique physical fluid 
with the same effective velocity. When T e7 < T < T- n the electron-photon rate dominates 
against the Coulomb rate which is anyway larger than the ion-photon rate. Finally for 
T < T C7 the Coulomb rate is always dominant which means that the ion-electron fluid 
represents a unique entity characterized by a single velocity which is customarily referred 
to as the baryon velocity. The effective temperatures T e7 and T e i determine the hierarchies 
between the different rates and should not be confused with the kinetic temperatures of 
the electrons and of the ions which coincide approximately with the photon temperature 
T 7 ~ T c ~ T ; . For instance after matter-radiation equality (T c — T 7 )/T 7 ~ 0(H/T ei ) and 
(T; — T e )/T 7 ~ 0(H/T ci ) where H is the standard Hubble rate at the corresponding epoch. 

Depending on the range of temperatures the effective evolution equations for the vortic- 
ities will change. In the regime T > T i7 the Coulomb rate can be neglected in comparison 
with the Thomson rates and the vorticities of photons, electrons and ions approximately 
coincide. For T e7 < T < T i7 the Ohm law can be easily obtained from Eq. (12. 5 p and it is 
given by 

£ + )J c x5 = - + -^-o 2 r 7e (£-t; e ), (2.26) 
o 3 pb e 

where it has been used that the baryon density pb = [vn\ + m e )n coincides approximately 
with the ion density in the globally neutral case and that n = a 3 no; furthermore, in Eq. 
(I2.26p . a denotes the electric conductivity 



to 2 I Aire 2 n e 

Anal ei V rn P a 



expressed in terms of the Coulomb rate and in terms of the electron plasma frequencj]^] 
By taking the curl of both sides of Eq. ( I2.26|) the following relation can be easily derived 



V x E + V x (v e x B) = + a 2 r 7e (aL - ou e ). (2.28) 

a 3 pb e 

Recalling now Eq. Q and (pDEj) . Eq. fl2T28|) becomes: 



dB -> _ -, V 2 B 4p 7 2 mi 
or Ana 3 pb e 



V x (v e x B) + - a 2 — r c7 (a; 7 - u e ). (2.29) 



5 The electron plasma frequency of Eq. (|2.27|) must not be confused with the vorticity 
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In the same regime the evolution equation for the vorticities of the ions and of the photons 
are, up to spatial gradients, 

PTi- — » 

d T Ui + Wdi = jd T B, (2.30) 

Pi a 4 

8 t uj 1 = ar 7e (u e — u; 7 ). (2.31) 

By eliminating the electron-photon rate between Eqs. (I2.30p and (I2.3ip and by neglecting 
the spatial gradients in Eq. ( I2.29p . the following pair of approximate conservation laws can 
be obtained 

d T ( auji + — B) = 0, (2.32) 



mi 
a 



d T [ — B - ) = 0, (2.33) 



where the ratio is given by: 



R h = = 30.36 (—)h 2 o n h0 . (2.34) 

4p 7 \ z J 

By further combining the relations of Eqs. f 1 2 . 3 2 j) and f !2.33|) the vorticity of the photons can 
be directly related to the vorticity of the ions since d T [R\,uj\ + cj 7 ] = 0. By assuming that at 
a given time r r the primordial value of the vorticity in the electron photon system is u r and 
that B(t t ) = we shall have that 

4 p (r r ) 

a r uji(r r ) H 7 r a r cJ 7 (r r ) = u v - (2.35) 

3pb(T"r) 

Thus the solution of Eqs. (12.321) and f !2.33|) with the initial condition f |2.35|) can be written 
as: 

— * / — * \ e B(x,t) a v _ 

Ui(x, t) = -— + — w r , (2.36) 

mi a{T) a{T) 

u; 7 (f, r) = Rh y^ [td T - a(r)uji(x, r)]. (2.37) 
a(r) 

The approximate conservation laws of Eqs. ( 12.32l) -( |2.33l) can also be phrased in terms of the 
physical vorticities Qx(x,t) = a(T)cu x (x,r) where X denotes a generic subscript^]. 

For typical temperatures T < T e7 the electrons and the ions are more strongly coupled 
than the electrons and the photons. This means that the effective evolution can be described 
in terms of the one-fluid magnetohydrodynamical (MHD in what folllows) equations where, 



9 Note that while ojx is related to B, the physical vorticity tlx is directly proportional to B. For instance, 
in the treatment of [21 EJ 0] the use of the physical vorticity and of the physical magnetic field is preferred. 



on top of the total current J the center of mass vorticity of the electron-ion system is 
introduced 

uiiUi + m e u e 

u h = ■ . (2.38) 

m e + m\ 

Equation (12 . 13[) (multiplied by m e ) and Eq. (12.151) (multiplied by m{) can therefore be 
summed up with the result that 



d T u h + Uu h = VX ( / X B) + - a h ). (2.39) 

a 4 p b 3 p b 



The evolution equation for the total current can be obtained from the difference of Eqs. (I2.5P 
and (12.61) . Since the interaction rates are typically much larger than the expansion rates the 
Ohm equation can be simplified and becomes 

— * — * J 4 777, ■ 

E + Vb x B = - + -^^a% e (v J -v b ), (2.40) 
a 6 pb e 

where v b is the baryon velocity related to the baryon vorticity as cJb = V x £f b . The similarity 
of Eqs. (I2.28P and f !2.39|) should not be misunderstood: while Eq. f !2.28|) follows from the 
right hand side of Eq. (12. 5ft . Eq. ( 12.391) follows by taking the difference of Eq. ( 12. 61) 
(multiplied by n() and of Eq. (12. 5p (multiplied by n e ). The expression obtained by means 
of the latter difference is rather lengthy and can be found in its full generality, in Ref. [13] 
(see, in particular, Eqs. (7) and (10)). Here the expression has been simplified by neglecting 
higher orders in (m e /m{). The effective set of evolution equations can then be written, in 
this regime, as 

V x (fx B) e! _ , 

d T u h + m h = \ J - + —(u^-u h ), (2.41) 

a Pb H\ 



b 



d T B = V x (v h xB) + + ^-e\Q h - w 7 ), (2.42) 

4-na e R\ } 

<9 T u; 7 = e'(uj h — w 7 ), (2.43) 

where e' = aT ei is the differential optical depth where, as usual, the contribution of the ions 
has been neglected. In the tight coupling limit Eqs. (I2.4ip . (I2.42p and (12.431) imply that 
uj hl ~ uj h ~ u-y while obeys 

Kith _ „ V x (fx B) 
d ^ + R^l^ = R % h a*(R h + iy (2 - 44) 

In analogy with what has been done before, the conservation laws can be derived by com- 
bining Eqs. (I2T4T|) and ( I2T42]) 



o/ - rrn _ \ -* ,_ V 2 B miV x (J x B) 

d T [B + — auj h )= Vx (v h xB) + - + — \ '-. (2.45) 

e / Aira e erpb 
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From Eqs. (I2.42p and ( I2.43P and by neglecting the spatial gradients it also follows 



d T [B - - u. 



'7 



) 



0. 



(2.46) 



Equations ( 12.451) and ( 12.46[) are separately valid, but, taken together and in the limit of 
tight baryon-photon coupling, they imply that the magnetic filed must be zero when the 
tight-coupling is exact (i.e. cJ 7 = u5b)- In spite of the various physical regimes encountered 
in the analysis of the evolution of the vorticity the key point is to find a suitable source of 
large-scale vorticity which could be converted, in some way into a large-scale magnetic field 
[28] (see also [29j|30]). The conversion can not only occur prior to matter-radiation equality 
but also after [5] in the regime where, as explained, the baryon-photon coupling becomes 
weak. Indeed, Eqs. (12.321) and (I2.45P have the same dynamical content when the spatial 
gradients are neglected and the only difference involves the coupling to the photons. 

There have been, through the years, suggestions involving primordial turbulence (see 
the interesting accounts of Refs. [21]), cosmic strings with small scale structure (see, e. g. 
[3"2~| |3"5| I3~4"]). Since matter flow in baryonic wakes is turbulent, velocity gradients will be 
induced in the flow by the small-scale wiggles of the string producing ultimately the vortic- 
ity. Dynamical friction between cosmic strings and matter may provide a further source of 
vorticity [33]. There have been also studies trying to generate large-scale magnetic fields in 
the context of superconducting cosmic strings (see, for instance, [M] and references therein). 
The possible generation of large-scale magnetic fields prior to hydrogen recombination has 
been discussed in [351 [361 EI] ( see a ls° [3H])- The vorticity required in order to produce the 
magnetic fields is generated, according to [35] , by the photon diffusion at second order in the 
temperature fluctuations. In a similar perspective Hogan [37] got less optimistic estimates 
which, according to [53 136] . should be attributed to different approximation schemes em- 
ployed in the analysis. Along this perspective various analyses discussed higher-order effects 
using the conventional perturbative expansion in the presence of the relativistic fluctuations 
of the geometry [39]. In the present paper, as already mentioned, we are going to follow 
a different route since we intend to use the gradient expansion for a direct estimate of the 
vorticity. 

3 Vorticity evolution in gradient expansion 

The conservation laws derived in section [2] hold under the hypothesis that the spatial gra- 
dients are neglected in the evolution equations of the vorticity. The logic of the gradient 
expansion [UJ [HI [191 I2Q1 [ZD [22] can be combined with the tenets of the drift approximation 
[23 1 [24 ] [25] m the context of the ADM decomposition [151 EE]- It will be shown hereunder that 
the resulting formalism [T4] provides a more general description of the angular momentum 
transfer between the various species of the plasma. Consider therefore the standard ADM 
decomposition where the shift vectors are set to zero but the lapse function kept arbitrary, 
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i.e. g 00 (x,r) = N 2 (x,t) and gy(x,r) = —jij(x,r). In this case the Maxwell equations can 
be written as 

d • E = 4ne[m - n e ], d ■ B = 0, (3.1) 



d T B + d x E = 0, d x B = Aire 



rii Vi - n e v e 



where the rescaled electric and magnetic fields are given by: 



+ d T E, (3.2) 



E\x, r) = (Xl) £<(x, r), B\x, r) = ( ^ ) B\x, r); (3.3) 

in Eq. (13.31) the subscripts specify that the rescaling is space-time dependent. The rescaled 
concentrations are 

ni(x,r) = y^yni(f,r), n e (x, r) = y/jn e (x, r). (3.4) 
The shorthand notation^ employed in Eqs. (I3.ip -( l3~3|) implies for a generic vector A 1 , 



8 ■ A = c^, (9 x Af = dj 



(3.5) 



In appendix El some relevant complements on this formalism have been collected to avoid a 
digression from the main line of arguments contained in the present section. Two relevant 
aspects must anyway be borne in mind: 

• in the conformally flat limit (i.e. N(x,r) — > o(r) and 7^(2?, r) — >■ a 2 (r)5jj) Eqs. (13. ip 
and ( 13. 2 p reproduce exactly Eqs. ( 12. 2 p and (\2.3\i : 

• the same comment holds for all the other fields (i.e. comoving or physical) involved in 
the fully inhomogeneous description. 

Using the generalized curl operator of Eq. ( 13 .5p the vorticity of the ions, of the electrons 
and of the photons can be written as 

where is the extrinsic curvature (see appendix |Aj while and A m are defined a 

A^ = NY k Y n Vnmk, AZ = 2N 2 [K ik y + 7 J 1 (3.7) 



10 Note that the operators introduced in Eqs. (I3.1I) - (|3.3I) are the generalized curl, divergence and gradient 
operators; they reduce to the conventional curl, divergence and gradient operators in the conformally flat 
limit. 

"Recall that rj abc = ^/je abc and that 7j abc = e abc /^j. 
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Using Eqs. (j3.6p -( )H77|) as well as the evolution equations of the velocities (see, Eqs. ( 1A.8[) — 
(1A.9j) ). the evolution for the vorticity of the electrons and of the ions can be written, respec- 
tively, ad3 



d T oj l c + [NK 
en P N 2 



N 



, i r>i , ,k 7~i 



d T u{ + ( A K 
en-.N 2 



4p 7 



(d x E) L + [d x (v c x B)f \ + NT^ut - ool) + -^JVT^K ~ w e), (3-8) 

" Pe 



, i ni , k 77! 



P. 



4p 7 



(d x E) % + [dx (tf, x E)f + ivr ie ^K - w*) + -^AT iT K - (3.9) 



Pi ' ' 3 pi 

Similarly, from Eq. ( 1A.10I) the evolution equation for the vorticity of the photons can be 
written as 



4 -NK - 
3 N 



cot - GiiZ -T= AT 7e (c^ - cot) + NT^ul - cot). (3.10) 



The quantities J 7 *, F{ and Tl. appearing in Eqs. (13. 8p . (13. 9p and (I3.10p are of the same order 
of the other terms appearing in the equations and they are defined as 



F l 



o l [X: i r;:>)+-Nr 7e d j 



Pj 

Pc 



+ d^Klvl - NdjKAZv? - d j 



fen e N 2 

V PeJl 



E m + (v e x B) 



(3.11) 



di ( A^ m ) + iNT+dj ( 9 ) A%K - v?) + Nd : 



4 



Pi 



Pi 



A lJ T- (v 



A' 



E m + (vi x £)' 



_ a N 2 ( A ij 

0, ( A!/ /■.*■)+ vftg* - -NdjK A k j v* - — d s \ -±d n 



P 7 7 



mfc 



(3.12) 
(3.13) 



7 ^9 3" K ~^ 1 

The generalized scalar and vector products appearing in Eqs. (13.111) . (13. 12[) and (13 . 13[) are 
defined as 

F ■ G = -f mn F m G n , (F x G) k = 7™^MF n G m r] ak , (3-14) 



N 

and coincide with the ordinary scalar and vector products in the conformally flat limit 
introduced after Eq. (13. 5p . The velocity fields appearing in Eqs. (I3.8P and (13 .9p are all 
subjected to the fully inhomogeneous form of the momentum constraint implying, from Eq. 

(HI), 

^ (ViK - V k K*) =£ 2 P (P + p)u°u i} u° = -^VTT^ 2 , (3.15) 



N 



N 



12 We shall focus, without loss of generality, on the situation where the lapse function is homogeneous, 
i.e. N(x,t) = N(t); in this case the already lengthy expressions will be more manageable since the spatial 
derivatives of the lapse function will vanish. 
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where u 2 = u % u^ij and where u° and u l can also be defined in terms of the total velocity field 
v l which turns out to be the weighted sum of the velocity fields of the electrically charged 
and of the electrically neutral species, i.e. 

(p + p)V k = E(Pa + P.K = PeVo + fWi + \p^ + -RuVl + p U*, (3.16) 
a 6 6 

where the contribution of the cold dark matter particles and of the massless neutrinos has 
been also added. The explicit connection between u°, u l and v l is given by: 

u = —— u=—coshy, cosh y = — (3.17) 
N N ^l-v 2 /N 2 

where v 2 = v l v^ij. In terms of v l and v 2 the momentum constraint of Eq. (13.151) can also 
be written as 

v 2 \ / 



%(p + p)jj=[l-jp)V k \K'»-K>y*). (3.18) 

All the discussion of section [2] can be generalized to the fully inhomogeneous case and 
we shall be particularly interested in the generalization of the conservation laws determining 
the angular momentum exchange between the various species. Consider then the situation 
where the electron-photon rate dominates against the Coulomb rate. In this case the fully 
inhomogeneous form of the Ohm law reads 



- E k - (v e x Bf + — + i- ^We, ^-{v* - v k e ) = 0. (3.19) 
a 3e pb N 



?k It* f>\k i i ^ H^, r V ' („,k „,fc 



By taking the generalized curl of Eq. ( 13.19)) (see Eq. (|3.5|) ) the following equation can be 
obtained 



-d x E - d x (v c x B) + 8 x ( J/tr) 

+ A ^m ; r e7 ^(a; 7 - a? e ) - | — ^(^ -« c )x9 
3e pb N 3 e 



r 



e7 



_7 P7 

iV p b 



0, (3.20) 



where, consistently with Eq. ( I3.14p . the last term at the left hand side is defined in terms 
of the generalized vector product and it vanishes exactly in the conformally flat limit. By 
assuming, as physically plausible prior to decoupling, that the conductivity is homogeneous, 
Eqs. ( 13. ip and ( 13. 2 p can be used inside Eq. ( I3.20p and the final equation will then be: 



d T B = d x (v c x B) — d x x B) — — —miT^ -^(w 7 — uj e ) 

Ana 3e pb N 

+\n 2 ^{v, - v e ) x d 



■p y 7 Pj_ 

67 N p b 



(3.21; 



Equation (13.211) reduces, in the conformally flat limit, to Eq. ( I2.29p . The same logic can 
be applied in all the other derivations and the obtained result expanded to first order in 
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the spatial gradients with the result that the generalized system for the evolution of the 
vorticities reads 

d T N" 



d T B k 



NK + 2- 



N 



ijj- 



4 r Pi 

■— r e7 — m 

3e 



3 AT 



1 jv ( 7 



en; 



c 7 



(3.22) 

(3.23) 
(3.24) 



Equations fl3T22|) . (13T23D and (13724]) reduce, respectively, to Eqs. (12T29D . ([2730]) and ([2731]) in 
the conformally flat limit. Equations (13 .22[) . ( I3.23[) and (I3.24[) apply in the situation where 
the magnetic fields are initially zero and do not contribute to the extrinsic curvature so that 
K\ = K/351+Kl with K\ = 0. In this case Eqs. ( I3.22[) ( T3.24[) reduce to a pair of remarkable 
conservation laws whose explicit expression, up to spatial gradients, is 



Or 



7 
N 2 

B k 



wf H -B k 

Pb 

mi 7 1 / 3 k 

=- - — U) 

eR h N 2 7 



0. 



0. 



(3.25) 
(3.26) 



where R^(x, Ti) is a constant in time (but not in space) and come from the inhomogeneous 
generalization of Rb(x, r): 

3p b (f,r) — 



Rb(x, t) 



i2 b (f,Ti)7 



1/6 



(3.27) 



4p 7 (f, r) 

The evolution of the vorticity of the baryons as well as the tight coupling between the baryons 
and the photons can be discussed in full analogy with the considerations already developed 
above in the case of the electron-photon coupling. The inhomogeneous generalization of the 
Ohm law when the Coulomb scattering dominates against both the electron-photon and the 
ion-photon coupling has been derived in Ref. jH] (see Eq. (3.34)). To leading order in the 
gradient expansion the evolution of the baryon vorticity can be written as 

d T N" 



d T U k : 

d T B k 



NK + 2 

771; e 



N 

7. 



(3.28) 
(3.29) 



e R h N 2 

where e' = iVT e7 is the inhomogeneous generalization of the optical depth. By eliminating 
e' between Eqs. ( 13.281) and ( 13. 29ft the following equation 



Or 



77 k , eni . 







(3.30) 



is readily obtained. Note that Eq. ( I3.30[) coincides, up to spatial gradients and in the 
conformally flat limit, with Eq. (12.45 j) . 
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4 Maximal vorticity induced by the geometry 

In this paper the the expansion is organized not in terms of the relative magnitude of the 
gravitational and electromagnetic fluctuations but in terms of the number of gradients carried 
by each order of the expansion. From the momentum constraint (see Eq. (I3.18P ). the total 
velocity field can be written, formally, 



S' 



NS i 
1 



1 - Vl + 4S 2 



NS* 



1-S 2 + 0{e 



+ 0{e 4 



£Up + p) 



ki 



(4.1) 
(4.2) 



where the orders of the expansion appearing in Eq. (14. ip are defined by the number of 
gradients. From Eqs. (14. ip and ( 13 . 6 [) — f j3TT|) the total vorticity can be written as 



tot 



dANAlS* 



1-S 2 + 0{e 



(4.3) 



To implement the gradient expansion let us parametrize the geometry as 

1 



Hi r) = a (r) (x) + tj (x, r)] , j lJ (x, r) 



o?(t) 



a v (x) -P v (x,t)}. (4.4) 



and keep the lapse function homogeneous, i.e. N(t) = a(r); otij(x) does not contain any 
spatial gradient while flij(x, r) contains at least one spatial gradient. The extrinsic curvature 
becomes: 



n 



K ik = — - 



H{a ik - 



a 2 a 

Furthermore we have also that the spatial Christoffel are: 



ka 



pm 
1 ab 



2 V a 



& ^nab + Ot A na b ~ P A 



nab 



-d T (3 ik 
2 ' 



(4.5) 

(4.6) 
(4.7) 



where X nab and X nab 



Xnab 9 n <y a i) ~\- 0f)Oi na -\- d a (Xb n , 

Xnab = —dnPab + d (3 na + Qafibn- 



(4.8) 
(4.9) 



The relevant term appearing in the momentum constraint becomes then 

a km d T d k /3 



V k [K km -K^ km 



V k K km + 



2a 3 



(4.10) 
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where 



V h K km 



+ 



2a 3 



i^{d a a)a am + Hd k a km + %——\ iab 
la I 

-d k d T (5 km - l -(d T (3 ab a mi y iab - (^)d T P° 



- Ud a Pa am + U 



d n Ot 



a 



(3 am + 2Ud k (3 km - Ua mi a ab \ iab 



+ H[a mi (3 ab + a ab p imi )A 



Equation (14. 10p can therefore be written as 



V k [K km - Kj km 



n 



^(d a a)a am + d k a km + ^ 
la 



mi „,ab 



-A 



iab 



jkm 



+ ^{a km d k d T (3 - d k d T (3 k 

'^d T f3 ab a™y tab +(^y T (3° 



+ u 



2d k f3 km - d a f3a am + 



d a a 



a 



ft 



am „ mi 



„ mi „ao \ 

a a X 



iab 



I „ ™ iQ a b i „ ab ami \ \ 

+ I a p + a p A 



iab 



The previous expression can also be recast in a more handy form: 



V fc (K 



r km 



km 



a 6 2a 6 



Z™(a,ft)-l™(a,ft)+Hl?(a,ft) 



where the three functionals of aiij(x) and 0ij(x,r) are defined as 

1 d a a 



Z m a) 



2 a 



-a am + d a a qm + 



a mq a ab 



qabi 



T^(a,/3) = a qm d q d T ft -d T d q ft qm } 

d a a 
2a 



X 2 m («,/3) 



a 



qui 



(d T ft ab )\ qab + ^d T (3 am , 

f) a 

, am i ^a^ nam 



7^(a, p) = 23 q ft qm - (d a ft)a am + ^f3 am + \ b ( a qm f3 ab + a ab ft mq 

a 



a 



mq _ ab \ 
OL X qab . 



(4.11) 



(4.12) 
(4.13) 

(4.14) 
(4.15) 
(4.16) 

(4.17) 



With the result of Eq. (14.131) we can compute the first relevant part of the final expression, 
namely: 



N^^VamnVklK^-K^"' 



a* 
a kj a m 



-Ha kl a iri Z m (a)t kr , 



+ T-L[a kj ft in + a in (3 kj )Z m (a)e kmn H — e k 



I?(a,f3)-Z™(a,ft) 



(4.18) 



16 



Recalling that, furthermore^! 

*>, s 9 37£?(1+U0 /ai\ 3u,+ 1 ,9 9 / \ 

tU P + p)a 2 = a{w+iy ^ l + p/ 2)w+1 (-j) , = 3W?. (4.19) 

Putting all the various parts of the calculation together we have that, from Eq. (14. 3p . 

u\ ot = djA ij , A ij = N2 ^+ V ? mn Va (K am - r m Kj , (4.20) 
then the quantity A^ becomes: 

(w+2)/2 , x 3w+l ( 

3nf(w + 1) Vai/ I 



+ ft 



a kj (3 in + a in /3 kj 



a kj a in 
Z '{ajekmn H — ^ 



- 2T(a,/3) + ftX 3 m (a,/3) 



Xr(«,/3) 

- j{w + 2)(3a k ^a in Z m (a)e kmn y (4.21) 



The first line at the right hand side of Eq. (14.21 j) does not contain any spatial gradient and 
it is therefore 0(a). The remaining part of the expression at the right hand side of the 
relation reported in Eq. (14.2ip are instead O{0). Sticking to the situation treated in the 
present paper the explicit form of (x, r) can be determined in terms of (x) by solving 
the remaining Einstein equations written in terms of the ADM decomposition [151 US]- For 
this purpose Eqs. (1A.5|) and ( 1A.7j) can be written, respectively, as 

Nf 2 

d T K - NTrK 2 = — ^ (3p + p) , (4.22) 

Nf 2 

d T K\ - NKKl - Nrf = -^(p - p)5f. (4.23) 
Using Eqs. (j4.4p and ( 14. 5 p into Eqs. (I4.22p . the following pair of conditions are obtained 

Sr-H = -^(/" + 3p m ). (4.25) 

To obtain Eqs. ( 14. 24ft and ( 14.251) the total pressure and the total energy density have been 
separated as: 

p(x,r) = P {0) (t) +p (1) (f,r), p(x,r) =p (0) (r) + p (1) (f,r), (4.26) 
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We shall assume that w, the dominant barotropic index of the fluid sources, is constant. 
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where p^(x, r) and p^(x, r) vanish in the conformally flat limit. Using Eqs. (14. 4p and (14. 5 p 
into Eqs. (I4.23[) two further equations are obtained and they are: 



^(^r) + nd t bl + ir^' + ari = ~^ p{1) ~ p(1))6 - (4 - 27) 

d T n + m 2 = -^( p (°) - p <°)). (4.28) 

Solving Eqs. (I4.25P and (I4.28P under the hypothesis of constant barotropic index (already 
assumed in Eq. (14.191) ). p^ and can be eliminated between Eqs. (14.241) and ( 14.2TI) and 
it turns out that 0ij(x,r) obeys the following evolution equation: 

d?jt + mdM + 51 (^^ T 2 /3 + 2^-Hd T ^ + 2a 2 rf = 0. (4.29) 

By solving Eq. ( I4.29[) the explicit form of {3y can be written in a separable form as {3j (x, r) = 
g(r)pl{x) where: 



g(r) = c? w+ \ (4.30) 

(4.31) 



H?(3w + 5)( y 3w + 1) 



' j,-* 3w 2 — 6w — 5 .,• 
™ + 4( 9lu + 5) Pix)ii 



Note that -P/(x) = r- (x, T)a 2 (r) accounts for the intrinsic curvature computed from aij(x). 
In Eqs. (14 .221) and (14. 23ft the contribution of the velocity fields and of the magnetic fields has 
been neglected because they are subleading to 0((3). In the following two sections we will 
therefore present the full estimate of the vorticity to first-order in the gradient expansion. If 
needed the first-order result, together with Eqs. (I4.30p and (I4.3ip can be used to estimate 
the vorticity to higher order. 



5 Vorticity to first-order in the gradient expansion 



The simplest parametrization of 
written as 



x) which does not contain spatial gradients can be 



-2* 



a 



det otij 



-6*(£) 



(5.1) 



In this case it is easy to show that Z m (a) = and therefore the first-order in the gradient 
expansion vanishes identically. In the ACDM scenario the scalar mode appearing in Eq. 
( 15. ip leads to a <C 1 and therefore, in practice, c%(if) is accurately estimated by 

5ij — 2^(x)5ij. To have a Z m (a) ^ the contribution of the tensor modes must be included 
will then given by: 



and OLij[x 



a tj {x) 



Sij + hij(x) 



a lJ (x) 



5 lJ - h tJ + h ik h{ 



a 



i - ~ A K h\ 



(5.2) 
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where is divergenceless and traceless, i.e. <9j/i lJ = h\ — 0. It must be borne in mind 
that the scalar and the tensor modes, in the ACDM scenario and in its tensor extension, are 
defined in terms of the conventional perturbative expansion. As a consequence of the latter 
statement, the informations on the spatial inhomogeneities of the model are not specified by 
assigning the analog oiij{x) (or 7y(a?,r) to a given order in the spatial gradients). On the 
contrary, as it is more natural, the scalar and tensor modes of the geometry are specified by 
assigning the corresponding power spectra at a given pivot scale. To evaluate the appropriate 
correlators defining the vorticity we shall need first to obtain the fluctuations in real space 
(as opposed to Fourier space). Therefore, as we will show in the present and in the following 
section, the idea will be first to compute the fluctuations in real space and then to use the 
obtained result for the determination of the correlators defining the vorticity. This procedure 
will circumvent the calculation of complicated convolutions and will also be perfectly suitable 
for the applications described in section |6j Using then Eq. (I4.14p we have that 

Z m {a) = d q a qm + a mq a ab d b a qa = h qm h ab d b h qa + h ap h b p d b h^. (5.3) 

From Eq. (I4.2ip the tensor A l:j (a, (3) can be computed to lowest order (i.e. by setting = 0) 
and the result will therefore be written, using Eq. (15. 3p . as 



V / a \ 3w+1 
Av ( a ) = [ _ ) e m v 



h a %d b h am + h mq h ba d b h q a 
Finally, the total vorticity can be derived directly from Eq. ( 14.20p 



+ 0{e 2 ). (5.4) 



-C(r,w)e mtJ d j 



h h^d b h am + h mq h d b h qa 



+ 0{e 



3\ 



T-L / a \ 3w+1 

£M = o * , J -) . (5.5) 



To give an explicit estimate of the primordial vorticity the relevant cosmological parameters 
will be taken to be the ones determined on the basis of the WMAP 7yr data alone [26l [27] . 
In the ACDM paradigm the sole source of curvature inhomogeneities is represented by the 
standard adiabatic mode whose associated power spectrum is assigned at the comoving pivot 
scale k p = 0.002 Mpc -1 with characteristic amplitude An 

-> 9-7T 2 -> / k \ n s — l 

(n(k,r)n(p,r)) = —V 1l (k)5^(k + p), V n (k) = A n [ v ) , (5.6) 

where n s denotes the spectral index associated with the fluctuations of the spatial curvature. 
According to the WMAP 7yr data alone analyzed in the light of the ACDM paradigm and 
without tensors modes [26l E7j the determinations of An and of n s lead, respectively, to 
An = (2.43 ± 0.11) x 10" 10 and to n s = 0.963 ± 0.014. The standard ACDM scenario, 
sometimes dubbed vanilla ACDM is defined by six pivotal parameters whose specific values 
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are, in the absence of tensor modes 14 ! 

(fi b , tt c ,tt de , h , n s , e re ) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088). (5.7) 

To estimate the correlation functions associated with Eqs. (j5.4p and ( 15. 5 p it is mandatory 
to know in detail the numerical value of the correlation function of the tensor modes of the 
geometry which have not been detected so far but whose specific upper limits will determine 
the maximal magnetic field obtainable from the vorticity of the geometry. The tensor modes 
of the geometry are described in terms of a rotationally and parity invariant two-point 
function 

(hijix, r) hiM r)> = J yP T (fc, r) (5.8) 

where the tensor power spectrum at the generic time r is given by the product of the 
appropriate transfer function multiplied by the primordial spectrum: 

V T (k,T)=M(k,k eq ,r)V T (k), V T (k) = A T ^-j ; (5.9) 

note that At is the amplitude of the tensor power spectrum and tit is the tensor spectral 
index. The transfer function Ai(k, k eq , r) can be computed under several approximations 
depending upon the required accuracy. The transfer function for the amplitude of the tensor 
modes can be numerically computed by solving the evolution of the tensor fluctuations across 
the matter-radiation equality and the result is [40j |JT] 



M(k, fc eq , r) ^ 9]!{kT) 



k \ ( k 



1+Ci[ ir +c \k 

rv e q / \ rv e q 



2i 



|Jbf " ' ~ M '~ 1 ' ~ M ' ' (5 ' 10) 

whereEl according to [401 [4Tj . c\ = 1.26 and c 2 = 2.68. In Eq. ( 15. 10}) jx(y) = (smy/y 2 - 
cos y/y) is the spherical Bessel function of first kind which is related to the approximate 
solution of the evolution equations for the tensor mode functions whenever the solutions are 
computed deep in the matter-dominated phase (i.e. a(r) ~ r 2 ). Instead of working directly 
with At it is often preferred to introduce the quantity customarily called r T denoting the 
ratio between the tensor and the scalar amplitude at the pivot scale k p 

= £ = (5.11) 
An Vn{k P ) 



14 Following the standard notations (slightly modified to avoid possible clashes with previously dehned 
variables) fib, ^c^dc denote, respectively, the present critical fractions of thebaryons, of the dark matter, 
of the dark energy; ho is the Hubble constant in units of 100km/(secMpc), n s is the scalar spectral index 
while e rc denotes the optical depth at recombination. 

15 The analysis of [32] gave c\ = 1.34 and C2 = 2.50 which is fully compatible with the results of [101 HTj . 
In the approach of [42] (see also [43]) the calculation of the amplitude transfer function, in fact, involve a 
delicate matching on the phases of the tensor mode functions. Conversely, if the transfer function is computed 
directly for the spectral energy density, the oscillatory contributions are suppressed as the wavelengths get 
shorter than the Hubble radius (see below). 
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In principle can be taken to be independent of and this possibility will also be con- 
templated in the present discussion. At the same time, if the scalar and the tensor modes 
are both of inflationary origin, ut is related to rx and to the slow-roll parameter e which 
measure the rate of decrease of the Hubble parameter during the conventional inflationary 
stage of expansion: 



Ut = — - 



r T 



-2e, 



H 



(5.12) 



the overdot denotes the usual derivative with respect to the cosmic time coordinate; in Eq. 
(I5.12p the spectral index is frequency-independent but there exist situations where more 
general possibilities can be contemplated such as, for instance 



ht = — 2e + — In {k/h 



1) + 



8 



(5.13) 



If «t = the tensor spectral index does not depend upon the frequency and this is 
the case which is, somehow, endorsed when introducing gravitational waves in the minimal 
tensor extension of the ACDM. If a tensor component is allowed in the analysis of the WMAP 
7yr data alone the relevant cosmological parameters are determined to be 



(fi b , tt c ,tt de , h , n s , e re ) = (0.0430, 0.200, 0.757, 0.735, 0.982, 0.091). 



(5.14) 



In the case of Eq. (15. 7p the amplitude of the scalar modes is An = (2.43 ± 0.11) x 10 9 
while in the case of Eq. (15.141) the corresponding values of An and of are given by 



A n = (2.28 ±0.15) x 10" 



r T < 0.36, 



(5.15) 



to 95 % confidence level. To avoid confusions it is appropriate to spend a word of care on the 
figures implied by Eqs. (I5.14p and (I5.15P which have been used in the numeric analysis just 
for sake of accuracy. The qualitative features of the effects discussed here do not change if, 
for instance, one would endorse the parameters drawn from the comparison of the minimal 
tensor extension of the ACDM with the WMAP 5yr data release [441 145] . implying, for 



instance, An = 2.1±|;§ x 10 



0.984 and rx < 0.65 (95 % confidence level). Similar 



orders of magnitude can be also obtained from even older releases [461 El- 



6 Magnetic field induced by the total vorticity 

The total vorticity derived in the previous sections is larger than the vorticity of the ions. 
Therefore, the total magnetic field derived on the basis of u\ ot is larger than the one derived 
on the basis of the ion contribution. Of course this statement holds in an averaged sense 
since what matters is not the vorticity itself but rather its two-point function which will be 
explicitly computed in the present section. Using Eq. ( 15.51) the maximal obtainable magnetic 
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field will be the one given by Eqs. f l3.25p - fl3.26p or f )3.30p where the total vorticity induced 
by the geometry is given by Eq. (15.51) 



BLJ$,t) = -M<4(^r). (6.1) 
e N A ri\ 



which can also be written, by explicitly keeping track of the number of gradients, as 



Bi^ix, r) = |£(r, w) e^dj ^h al h\d b h am + h mq h ba d b h^ + 0(e 3 )} a(r 



1 + 0{e 2 



• (6-2) 



The prefactor appearing in Eq. (16.11) has been estimated, in Eq. (16.21) . by recalling that, to 
lowest order in the gradient expansion 

d T pi = NK Pi , d T Ui = NKh u (6.3) 

implying that p ; and hi scale in the same way with ^7 since NK = —d T In ^fiy. But then, 
from Eqs. flOO]) . f TOT]) and : 



1 + 0(e 2 ) 



(6.4) 



where the first correction is 0(e 2 ) and depends on (3 (see, e.g. Eqs. (I4.30p and (I4.3ip ) but 
it will be immaterial for the present ends. From now on the subscripts will be dropped 
but it will be always understood that we are referring here to the total vorticity and to the 
maximally achievable magnetic field. As a consequence of Eq. (16. ip the correlation function 
of the magnetic field can be related to the correlation function of the vorticity. To estimate 
the correlation of the vorticity and to obtain an explicit expression the key point is to reduce 
the six-point function of the tensor modes to the product of two-point functions. For this 
purpose it is not sufficient to consider the trace of the two-point function introduced in Eq. 
(I5.8P but it is rather necessary to proceed with the full tensorial structure of the correlator 
whose general parity and rot at ionally- invariant form will be denoted as 

(x,t) h mn (y,r)), (6.5) 

where Gij mn (r) is only function of r = \r\ where r = x — y. Since both hij(x, r) and h mn (y, r) 
are transverse and traceless, Gij mn (r) will have to share the same properties. In particular, 
Gijmn{f) must be symmetric for i — > j, m — > n, (ij) — > (mn) and satisfy the following 
properties 



d_ 

dr 

dk„ ,,. sin kr 



Gijmn 0; G ' %% m n Gij rnm (6.6) 



r- ^ -i j f ClrC / ■, \ Sill rCV , „ 

Tr[G ijrm ] = Gijij = J — "Pt(^) —j^r- (6.7) 
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The properties of Eq. ( 16. 6 p and ( 16. 7p are a reflection of the divergenceless and traceless 
nature of hij(x, r) while the requirement on the trace follows from the consistency with Eq. 
(15. 8p . The general form of Gij mn can therefore be written as 

Gijmn( r ) — {^im^nj + ^mj^ni^j G\(r) + 5ij5 mn G^fV) 

+ (Sij r m r n + 5 mn Ti r^jG z {r) 

^^jn^i^m ~\~ &vmTjT n ~\~ 5j m TiT n -\- 5i n T jTfy^jG 4(?") 

+ rirjr m r n G 5 (r), (6.8) 

where the various independent functions appearing in Eq. (16. 8p are determined in appendix 
IB1 The methods used to analyze the real-space correlators are the ones exploited in usual 
applications of statistical fluid mechanics [HI 09] . To evaluate Eq. (16.21) can then proceed 
as follows. Using Eq. (15. 5p . the explicit form of the correlator of the vorticity becomes 



Thmh' m.' 



h a ' V h a > q>h q f m r 



). (6.9) 



ha b h a q hq m 

III L Ay 

By defining (wHr)) = (co l (x, r) u l (y, r)) and by recalling the notations of appendix |B| we 
shall have thatl 16 ! 

<u, 2 (r)> = C\r, W )e^e^^^^T M (6.10) 

where the quantity T m >b> m b is a function of r; the explicit form of T m 'b'mb is given in appendix 
[B]in terms of the two-point functions Gij mn . Furthermore, since Tm'b'mb can be written, in 
general terms, as 

Tm'b'mb = Tx^^mibSmb' + 5 mm '5bb') + T2(r)5 m 'b'S m b 

+ T 3 (r)(5 m > b 'r m r b + 5 mb r m >r b >) + T A (r)(Sbm'r m rb' + 5 mb'^b^m' 

+ Smm'nrv + 8wr m r m >) +r m >rvr m r b T 5 (r). (6.11) 

By using the results of appendix [B] the explicit values of the five T»(r) can be expressed in 
terms of the two-point functions Gij mn (see Eq. (IB. 22ft and (1B.231) (1B.27I) ). There are two 
physically complementary regimes where the primordial vorticity and hence the magnetic 
field can be evaluated. Comoving lengths tq defined between 1 and 100 Mpc are smaller 
than the Hubble radius at equality sine 



l "MO 



16 Recall that r = x — y, i.e. r = \r\ = \o£— y\- 

17 The quantity rx (denoting, in sec. [5j the tensor to scalar ratio) must not be confused with tq and r eq 
which denote specific values of the radial coordinate. 
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where H is the present value of the Hubble rate, f2 M0 is the present value of the critical 
fraction in matter and Quo is the present value of the critical fraction in radiation. The pivot 
length r p = 500 Mpc at which the tensor amplitudes are assigned is such that r G < r eq < 
r p . Therefore, after matter-radiation equality and, in particular, at photon decoupling, the 
correlation function of the magnetic field can be estimated as 

(E 2 (r)> = 6.348 xloW-^-Vf ^_\ 3 f z ^ + 1 Y 

v v ;/ V0.32/ V2.43 x 1Q- 9 



1089.2 

C(n T ,r) G 2 , 



0.134 J V4.15 x 10- 

C(n T ,r) = c(n T ) - +d(n T ), (6.13) 



in units of G 2 = Gauss 2 and where the constants c(nx) and din^) are given byf^l 
c(n T ) = -2(n T -4)(n T -3)[2r2 T (r2 T -6) + 19]cos 3 (^^jr 3 (n T -5), 
i(nT) = 36+14^-4) 

The typical values of are negative and C(10~ 2 ). Indeed, assume, consistently with Eq. 
(I5.15p . that r T ~ 0.32. Then, according to Eq. (I5.12p . n T ~ —0.04 and e ~ 0.02. Concerning 
the results of Eqs. (16.1 3D and (I6.14p few comments are in order: 

• the prefactor C(t,w) is estimated in the hypothesis w = 0, a\ — a eq and "Hi = H eq 
since we ought to estimate the field prior to photon decoupling; 

• recalling that 7-L = aH the value of the Hubble rate at the equality time can be 
estimated as: 

H eq = y^tfo (^) 3/2 = 1-65 x 10~ 56 (f^) 2 M P ; (6.15) 



• the result of Eq. (I6.13P holds for comoving scales r < rp = 500 Mpc (which are the 
ones relevant for the gravitational collapse of the protogalaxy) and it is not sensitive 
to the variation of r provided np is nearly scale-invariant; 

• if r T ~ 0.32, then n T = -0.04; from Eq. (JSZHD, C(r G , -0.04) ~ 0.07 while for 
r = 100 r G we have that C(100r G , -0.04) ~ 0.01. 

By thus approximating C(nx,r) ~ 0(1) in the range r = 1-100 Mpc and for 0.2 < < 0.3 
we get the following value for B max = J (B 2 (r)) 

r T \ 3 /V M \^ 2 fz dec + Y 



B - 2 519xl0- 38 r T V7 An X ( 
^ max " 2 - 5i9xiU U.32J l2.43xlO^J I 



1089.2 



ft °" R ° ,)"*G. (6.16) 
0.134 J V4.15 x 10-V V ; 



18 Recall that because of the relation (|5.12[) tit < and rT > 0. 
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The result of Eq. (16.161) does not seem to be even remotely relevant for galactic magnetoge- 
nesis or for cluster magnetogenesis. In spite of the intricacy and of the ramification of the 
galactic dynamo hypothesis, it is useful to compare Eq. (I6.16P with the minimal requirements 
stemming from what we would call optimal or ideal dynamo, namely a process where the 
kinetic energy of the protrogalaxy is converted into magnetic energy with maximal efficiency 
Let us denote with N mt the number of (effective) rotations performed by the galaxy since 
gravitational collapse and with p a and pb the matter density after and before gravitational 
collapse . 

The typical rotation period of a spiral galaxy is of the order of 3 x 10 8 yrs which should 
be compared with 10 10 yrs, i.e. the approximate age of the galaxy. The maximal number of 
rotations performed by the galaxy since its origin is then of the order of N TOt ~ 30. Under 
the hypothesis that the kinetic energy of the plasma is transferred to the magnetic energy 
with maximal efficiency, the protogalactic field will be amplified by one efold during each 
rotation. The effective number of efolds is however always smaller than 30 for various reasons. 
Typically it can happen that the dynamo quenches prematurely because some of the higher 
wavenumbers of the magnetic field become critical (i.e. comparable with the kinetic energy 
of the plasma) before the smaller ones. Other sources of quenching have been recently 
discussed in the literature (see, for an introduction to this topic, section 4.2 of [50J and 
references therein). There is also another source of amplification of the primordial magnetic 
field and it has to do with compressional amplification. At the time of the gravitational 
collapse of the protogalaxy the conductivity of the plasma was sufficiently high to justify the 
neglect of nonlinear corrections in the equations expressing the conservation of the magnetic 
flux and of the magnetic helicity. The conservation of the magnetic flux implies that, during 
the gravitational collapse, the magnetic field should undergo compressional amplification, 
i.e. the same kind of mechanism which is believed to be the source of the large magnetic 
fields of the pulsars. Taking into account the two previous observations the estimate of Eq. 
( I6.16P must be compared with the bound 



in nG units. Even assuming N TOt = 30, p a ~ 10 24 g/cm 3 , and pb — 10 29 g/cm 3 the 
minimal value of -Bbound is C(10~ 25 )G. Clearly, by comparing Eq. ( 16.161) with Eq. ( 16. 17ft . 



Going then to cluster magnetogenesis, the typical scale of the gravitational collapse of a 
cluster is larger (roughly by one order of magnitude) than the scale of gravitational collapse 
of the protogalaxy. The mean mass density within the Abell radius ( ~ 1.5/ig" 1 Mpc) is 
roughly 10 3 times larger than the critical density since clusters are formed from peaks in the 
density field. Moreover, clusters rotate much less than galaxies even if it is somehow hard to 
disentangle, observationally, the global (coherent) rotation of the cluster from the rotation 
curves of the constituent galaxies. By assuming, for instance, N rot = 5, a density gradient 
of 10 3 and 500 nG as final field, Eq. ( 16.171) demands and initial seed of the order 0.15 nG. 




(6.17) 



B, 



max 



"C -Bbound- 
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Another application of the results obtained in the previous sections can be the estimate of 
the magnetic field induced by the total vorticity for scales which are larger than the Hubble 
radius prior to matter radiation equality. To conduct this estimate the explicit form of the 
correlators will change. First of all in the pre-factor £(t,w) we shall choose w = 1/3 and 
ai = a r and %i = 7-L T with H T ~ 10~ 5 M P . Thus for typical length-scales larger than the 
Hubble radius at equality and for typical times of the order of the equality time the analog 
of Eq. (I6.13P can be written as 

/ f \ ~4-3n T 

C(n T ,r) = c{n T ) - + d(n T ), (6.18) 

\Tp/ 

where the numerical constants c(tit) and d(nx) are given by 

c(n T ) = -2n T (n T + l)[2n T (n T + 2) + 3]cos 3 ^^^r 3 (n T -l), 

= _ 2(7,| + 28nx + 18) 
v ; 45^(n T + 2) 2 (n T + 4) v ' 

Equation (16.181) holds under the assumption r < r p which means, in practice, that it applies 
only for a narrow range of scales 120 Mpc < r < 500 Mpc. If r ~ 250 Mpc then r/r p = 0.5 
and C(nT,r) ~ (9(163) and 

V0.32/ V2.43X10" 9 / V 0.134 J K 1 



7 Concluding remarks 

The idea explored in this paper has been to compute the vorticity by employing a recently 
devised framework for the treatment of fully inhomogeneous plasmas which are also grav- 
itating. The latter description brings a new perspective to the study of the evolution of 
the vorticity exchange in the electron-ion-photon system without postulating the customary 
separation between a (preferably conformally flat) background geometry and its relativis- 
tic fluctuations. A set of general conservation laws has been derived on the basis of the 
fully inhomogeneous equations in different temperature regimes depending on the hierar- 
chies between the exchange rate of the vorticity between electrons, ions and photons. After 
expanding the Einstein equations as well as the vorticity equations to a given order in the 
spatial gradients, the total vorticity has then been estimated to lowest order in the gradient 
expansion. 
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The maximal comoving magnetic field induced in the ACDM paradigm depends upon 
the tensor to scalar ratio and it is, at most, of the order of 10 G over the typical comov- 
ing scales ranging between 1 and 10 Mpc. The obtained results are irrelevant for seeding a 
reasonable galactic dynamo action and they demonstrate how the proposed fully inhomo- 
geneous treatment can be used for a systematic scrutiny of pre-decoupling plasmas beyond 
the conventional perturbative expansions. The estimate of the primordial vorticity induced 
in the ACDM scenario can also turn out to be relevant in related contexts such as the ones 
contemplated by non conventional paradigms of galaxy formation. 
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A Gradient expansion and pre-decoupling physics 



In this appendix we are going to recap the essentials of the fully inhomogeneous description 
of pre-decoupling plasmas already introduced in Eqs. (I3.ip -( l3~5"l) . We will follow here the 
formalism developed in Ref. [2] and describe the fully inhomogeneous geometry in terms 
of the ADM decomposition [151 116] : 
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N k N k , 



9ij = 
N*Ni 
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ij 
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In the ADM variables the extrinsic curvature Kij and the spatial components of the Ricci 



tensor become: 
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2iV 

a (3)-pm a(3)-pm _i_(3) pm (3)-pr. 

u m 1 in U n 1 im. ' 1 ij 1 r 



-drla + (3) ViNj +< 3 > 



(3) pm (3)pn 



J ii 



(A.2) 
(A.3) 



Defining the total energy-momentum tensor of the fluid sources, the contracted 

form of the Einstein equations reads 



T 

rpu 



T = ^T. 



rpfl 



(A.4) 



As in the bulk of the paper we are now going to focus on the situation where the shift vectors 
vanish and the lapse function is homogeneous but time dependent (i.e. N(x,r) = N(t)). 
The (0 0), (ij) and (Oi) components of Eq. (1A.4|) become then: 



d T K - NTiK 2 + V 2 7V = N%>&±£ + (p + p) u 2 }, 



ViK - V k K k = N£ 2 P u° Ui(p + p), 
d T K( - NKK( - Nr{ + VN j N 



P-P 



5\ - (p + p)uiU 3 ], 



(A.5) 
(A.6) 
(A.7) 



where u 2 = u % v? r yij. The electron and ion velocities appearing in Eqs. (13. ip and (13. 2 p reduce, 
in the conformally flat case (i.e. N(t) —> a(r) and 7y(af, r) — > a 2 (r)5ij) to the velocity fields 
appearing in Eqs. ( 12. 3ft . ( 12. 13ft and (12. 15ft . In the fully inhomogeneous case, the evolution 
equations for the velocities of the electrons, ions and photons can be written, respectively, 
as 



d T v k p + Nd k N - g k vi 



eh e N 2 



'3 "e 



PcV7 

k 



d T v k + Nd k N - Q k v{ 



+ NT ei (v* 
ehiN 2 



E k + [y e x B) k 
.M , 4 Pi 



3 Pc NT ^ 



(A.8) 
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d T v* + Nd k N - 



+ jv^r to ( v ;-«*) + ^Jvr 

Pi 3 pi 



(A.9) 



4p 7 ' 

+ ivr 7e (^ - v k ) + ivr 7i (^ - «*), 



(A.10) 



where 



(A.ll) 



As in the conformally flat case the evolution equations of the electrons and of the ions 
can be combined by defining the center of mass velocity of the electron-ion system = 
(m e v k + rriiV k ) / (m e + rrii) so that the the effective evolution equations for the baryon-lepton- 
photon fluid become 



<9 r p 7 

drV* 



fh k 
N li- 



ck NK zk 



7p b (l + m e /m i ) 3p b ' 
N 2 / \ 
V ^ ~ 1i Qm Y^) ~ N9kN + NT ^ ~ 



(A.12) 
(A.13) 
(A.14) 



where and p 7 denote, respectively, the photon velocity and the photon energy density. 



B Some relevant correlators 



The correlator appearing in Eq. (16.91) . i.e. 



bmb'm' 



hab haq hqm 



h a ib'h a i q'hq'm 1 



(B.l) 



must be computed in terms of the corresponding two-point functions in real space (see Eq. 
(16. 5p ). The general form of the two-point function in real space has been already mentioned 
in Eq. (16. 8p and the functions Gj(r) (with % = 1, ... 5) are given by: 
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(B.3) 
(B.4) 
(B.5) 
(B.6) 
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where Fi(r), F 2 (r) and F 3 (r) are fully determined once the power spectrum is known and 
are defined as: 



, i 1 f dk _ , , . sin kr „ , , 1 f dk ^ 
dk 
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Using Eq. (lR7l) inside Eqs. fjR2]) - (jET6]) we have that 
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where jo(kr) and ji(kr) are spherical Bessel functions of zeroth and first-order [5TJ 
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It is useful to compare the two different asymptotic limits of the various Gi(r), i.e. for kr < 1 
and for kr > 1. In the limit fcr > 1 we have that: 
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Conversely, in the limit kr < 1, Eqs. ( 1B.8f) -( lB.12l) imply: 
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The explicit form of the two-point function Gij mn implies that the six-point function appear- 
ing in the correlator of the vorticity can be expressed as 

5 

Tbmb'm' = ^ Tbmb'm'i (B.22) 



where the 5 distinct contributions correspond to 
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The overline signifies that the corresponding correlator is evaluated in the limit r — > 0. 
According to Eqs. (!B~T7l) -( lB~2ll this limit is non-singular. Notice finally that in terms of 
Tbmb'm' the correlation function of the magnetic field can also be written, with shorthand 



notation, as 
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The real space approach is more effective and convenient for an explicit estimate of the 
vorticity and the idea is therefore to express the correlation functions in real space, take 
the appropriate derivatives and then expand the result in the desired limit. Denoting with 
R = r/r p and with x = k/k p , the integrals over k appearing in Eqs. (1B.8l) -( lB.12l) can be 
computed explicitly by changing variable and by using the following pair of relations [53] : 

R 2 - 



x n m sin xRdx 



x 



cos xR dx 



R 

m — n — 2 
+ R m ~ n - l cos 
1 



1F2 



ai]bi,b 2 ; — 



TC(m — n) 



r[l — m + n] 



m — n — 2 



1F2 



a 1 ;b 1 ,b 2 ] — — 



+ R 



m—n—l 



sin 



7r (m — n) 



r[l — m + n], 



[B.29) 



(B.30) 



31 



where n < m (i.e. (n — m) is negative). In Eqs. (IB.29[) and (IB.3Q[) p F q [ai....a p ; b±, ...b q ; z] 
denotes the generalized hypergeometric function of argument z; in the case of Eqs. ( IB. 291) 
and flB~30|) , p = 1, q = 2 and 
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The integrals are taken from 1 to infinity since the integral over k starts from k p implying 
that the lower limit of integration in x is 1. Equations (1B.29j) and ( ]B.30j) can be used to derive 



the real space form of the correlator of Eq. (16.51) . Using Eqs. ( 1B.29j) and (1B.30|) . the explicit 
form of Eqs. (1B.8I) - (1B.12|) can be derived and inserted into Eq. (16. 8p whose explicit form 
determines the real-space expression of the two-point functions of the vorticities through 
Eqs. ( 1B.23I) - (1B.27I) . After taking the appropriate derivatives the obtained result can be 
expanded in the wanted limits (e. g. R 3> 1 or R <C 1). The explicit real-space expressions 
of Eqs. (lB.8l) - f)B.12j) are typically rather lengthy but they are conceptually straightforward. 
This is why they will be omitted and only an example will be given. Even if the scales 
interesting for section [6] will be the ones close to the galactic scale, consider, for instance 
the expression of Gi(R) in the opposite limit, i.e. comoving scales much larger than r eq . In 



this case the expressions simplify since Ai(k, k eq , r) — > 1. Therefore, using Eqs. (1B.29j) and 
(lB~30l) . Eq. flEHD becomes: 
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To conclude this appendix let us show that the expression of Gij mn given in Eq. (16.51) 
coincides with the result directly obtainable in the particular case where the tensor modes 
of the geometry are quantized in terms of gravitons. When hij (x, r) is a field operator its 
expression can be written as [ffij 

\/2£p 



hij(x, r) 



(27r) 3 /2 a ( r ; 



E 

A 



d A k e%\k) 



h,x{r)e- ik - S + a\ x r k}X {r)e 



ik-x 



(B.34) 
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which also implies, using the properties of the creation and annihilation operators, 

dk 



G ijmn (r) = (hij(x, r) h mn (y, r)) = / —V T (k,T) Q ijmn (k) j (kr 



where 



Mk,r) 

Qijmn 



\fk{r)\' 



p 

AP 2 — — — 



^ ij mn a 



Pmi{k)P n j{k) + P m j [k)P n i[k) Pij (A;)P mn (fc) 



(B.35) 



(B.36) 
;(B.37) 



Pij(k) = (Sy — kikj) with k % = k l /\k\. In Eq. (1B.35j) it has been used that 

(0\a p ^al x \0) = 5^(k-p)5 x ,. 



(B.38) 



Furthermore, to derive Eq. (IB. 371) . the explicit form of the two tensor polarizations can be 
written, in explicit terms, as 



(B.39) 



where a, b and k are three mutually orthogonal unit vectors. 
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